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ABSTRACT 

In this paper, the variational iteration method was used to obtained numerical solution of Michaelis- 
Menten nonlinear biochemical reaction system. The correctional function for the system of equations were 
obtained and consequently the solution was obtained iteratively. Results obtained were compared with that of 
differential transform method (DTM) and were found to be in agreement. 
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I. INTRODUCTION 

The variational iteration method as proposed in [1,2] which is a modified general Lagrange multiplier method[3] 
has been shown to solve effectively, easily and accurately, a large class of nonlinear problems with approximation which 
converges quickly. It was successfully applied to autonomous ordinary differential equation [4], nonlinear partial differential 
equation with variable coefficient[5], Schrodinger-KDV, generalized KDV and shallow water equations[6], Burger's and 
coupled Burger's equation[7], the linear Helmotz partial differential equation [8] and recently to nonlinear fractional 
differential equation with Caputo differential derivative[9]. 

1.1 Mathematical model 

Consider the model 

E + A^Y^E + X (1) 
As provided in [10] 

Where E is the enzyme, A the substrate, Y intermediate complex and X the product. The time evolution of 
equation (1) can be determined from the solution of the system of coupled nonlinear ordinary differential 
equations 

dA 

= -K t EA + K ,y , (2) 

dt 

dE 

dt 
dY 

dt 
dX 



- K j£A + (K , + K 2 )Y , (3) 
K ^EA - (K l + K 2 )Y , (4) 
K 2 Y, (5) 



dt 

Subject to the initial conditions 

A(0) = A ,E(0) = E ,Y (0) = Y ,X(0) = X Q where the parameters K , , K _ 1 , A* 2 are positive rate constants 

for each reaction. System (2)-(5) can be reduced to only two equations for A and Y and in dimensionless form of 
concentrations of substrate, x, and intermediate complex between enzyme and substrate, y , [10] 

dx 

— = -x + (P - a)y + xy (6) 
dt 

dy 1 

— = -(x - py - xy) (7) 
dt e 

Subject to initial conditions 

x(0) = 1, y(0) = for the case a = 0.375 ,/? = 1.0 and e = 0.1 (8) 
Where a , ]3 and e are dimensionless parameters. 
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II. VARIATIONAL ITERATION METHOD 

Inokuti[3] proposed a general Lagrange multiplier method for solving nonlinear differential equations 
where he used nonlinear operator of the form as follows: 

lu + nu = g (x) (9) 
Where L is a linear operator, N is a nonlinear operator, s ( x ) is a known analytic function and U is an 
unknown that to be determined. 

The Inokuti method was modified by [4] which can be written as 

U n+l (x )= U „(*„) + p( L ^„ + N U n - g)ds (10) 



Where jj is an initial approximation and U n is a restricted variation. For arbitrary x 

U n + i {x)= U „(*)+ \*{LU n (s)+ NU n (s) - g(s))ds (11) 
o 

The above integral in equation (11) is called a correctional function and index fl denotes the nth 
approximation. Also equation (11) is called variational teration method (VIM). With Lagrange multiplier 
A [ii] 

(-1)> -x)" _I 

A = — (12) 

(m - l)! 

Where m is the highest order of the given equation. 
3 Application 

Consider equations (6) and (7). Using equation (11), the correct functional for equation (6) and (7) are 
X n+1 = X n + jUXl + X n - (P -a)Y n - X YJdx 

(13) 

y „+i = Y „ + \^ Y 1 - -( x „ - P Y „ - x JJ^y 

The Lagrange multiplier from (12) is 

A = -1 (14) 
Solving equation (13) and equation (14) we have 

X, = 1 - t 

(15) 

Y 1 = 10 t 

17 .25 f 2 10 t 3 

X = 1 - t + 

2 3 

(16) 

, 100 t 
Y 2 = 10 t - 105 t + 

3 

X (0 = 1 - t + 8.62500000 t 2 - 63 .08333333 t 3 + 62 . 1 8750000 r 4 - 194 .45833333 t 5 + 106 .25000000 r' 
- 15.87301587 f 7 

F(0 = 10/ - 105 r + 762 .08333333 r 3 - 819 .79166670 t" + 2077 .91666700 f 5 - 2020 .83333333 ? 6 
+ 476 .19047620 t 1 

(17) 

The numerical solution above are in agreement with that of the differential transformation method 
(DTM) stated by [10] which clearly indicate variational iteration method as a powerful tool. 
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III. CONCLUSION 

Variational iteration method was presented for Michaelis-Menten nonlinear biochemical reaction system. 
The method was found to be easy and efficient. Result obtained were compared with differential transformation 
method and are in agreement. The method is found to be accurate. 
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